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We elucidate the physical basis for the upper bound on high energy neutrino fluxes implied by the 
observed cosmic ray flux. We stress that the bound is valid for neutrinos produced either by p, 7 
reactions or hy p—p{n) reactions in sources which are optically thin for high energy protons to photo- 
meson and nucleon-meson interactions. We show that the upper bound is robust and conservative. 
The Waxman-Bahcall bound overestimates the most likely neutrino flux by a factor ~ 5/r, for small 
optical depths r. The upper limit cannot be plausibly evaded by invoking magnetic fields, optically 
thick AGNs, or large hidden fiuxes of extragalactic protons. We describe the implications of the 
bound for future experiments including the AMANDA, ANTARES, Auger, ICECUBE, NESTOR, 
and OWL/AIRWATCH detectors. 
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I. INTRODUCTION 

The observed cosmic ray flux at high energies implies 
an upper bound on the high-energy neutrino flux pro- 
duced in astronomical sources that are, like gamma ray 
bursts and the observed jets of active galactic nuclei, op- 
tically thin to photo- meson interactions on protons [y. 
The upper bound also applies to p — p or p — n collisions 
that create neutrinos via pion production in sources that 
are optically thin to nucleon-meson interactions. For sim- 
plicity, we shall refer in the following only to p ~ p colli- 
sions when we really mean either p — poip — n collisions. 

The high energy protons that produce neutrinos by 
p — 7 or by p — p interactions contribute to the observed 
cosmic ray flux after they leave the site where the neu- 
trinos are created. Since they have maintained a high 
energy in spite of the possibility of interactions with the 
cosmic microwave background, protons that are detected 
at earth with energies greater than 10^^ eV must orig- 
inate at redshifts z < 1. Therefore, the observed flux 
of high-energy cosmic rays determines the rate at which 
particles of those energies are being created in the rela- 
tively local universe. One can obtain an upper limit to 
the cosmic ray production rate in the whole universe by 
assuming that the rate increases with redshift like the 
fastest known population of astronomical sources, quasi- 
stellar sources (QSOs). Knowing an upper limit to the 
universal proton production rate, one can readily com- 
pute, by standard particle physics techniques, an upper 
limit to the rate of production of neutrinos by the same 
protons. 

The conservative upper bound derived in this way has 
come to be known in the literature as the Waxman- 
Bahcall upper bound ||l|. 

The upper bound exceeds what one can reasonably ex- 
pect to measure. As we shall see in Sec. Ill, the limit is 



derived numerically by assuming that all of the energy 



of the high energy protons produced in the astronomical 
sources is transferred to pions by photo-pion or nucleon- 
proton interactions. Each time that a 7r+ is produced by 
photo-pion interactions it receives only about 20% of the 
initial proton energy [nl . For p — p collisions at very high 
energies, the explicit calculation of the fraction of energy 
that is transferred to pions has not yet been done but is 
likely to be even less than for p — p interactions. Thus, at 
the very least, the Waxman-Bahcall upper bound is con- 
servative by a factor of five, and the most likely neutirno 
flux is about 5/t smaller than the Waxman-Bahcall limit 
for small optical depths. 

The Waxman-Bahcall bound is consistent with our pre- 
diction of the expected flux of high energy neutrinos from 
gamma-ray bursts (GRBs) Q,^. Naturally, the flux ex- 
pected from GRBs is less than the maximum allowed 
by the bound. For the GRB calculation, we estimated 
that ~ 20% of the total proton energy is transferred by 
photo-mesonic interactions to pions (proton-nucleon in- 
teractions are less efficient B ) . As described in the previ- 
ous paragraph, we assumed, in deriving our upper bound, 
that 100% of the proton energy is transferred to pions. 

The upper bound we set from the observed cosmic ray 
flux is two orders of magnitude lower than the intensity 
predicted in some previously published models for the 
production of neutrinos in AGN jets, implying that a 
km^ neutrino detector would record at most 1 neutrino 
from AGN jets per year. The same argument also rules 
out models in which most or all of the gamma-ray back- 
ground is produced by photo-meson interactions in AGN 
jets. 

Can the upper bound on the high energy neutrino flux 
be evaded? This is an important question since very 
large neutrino detectors are being designed for installa- 
tion in the ocean or a deep lake (see, e.g., Refs. Pfl), 
under Antarctic ice (see, e.g., Ref. 0), in space (see, 
e.g., Ref. [g|), and large area ground arrays (see, e.g., 



Refs. [| 10|). The design characteristics for these neu- 
trino detectors are being determined in part by the best 
available theoretical models. It is therefore necessary to 
examine carefully the justifications for different predic- 
tions of high energy neutrino fluxes. We do so in this 
paper. 

The goal of this paper is to clarify how the Waxman- 
Bahcall upper bound follows from the measured cosmic 
ray flux and, in the process of the clarification, to demon- 
strate that the upper bound is robust and conservative. 
We emphasize the implications of the bound for future 
neutrino and gamma-ray observations. 

There are two special types of sources for which the 
Waxman-Bahcall bound does not apply; these sources 
could in principle produce a neutrino flux exceeding the 
Waxman-Bahcall limit. The first special type of source is 
one in which neutrinos are produced by processes other 
than photo-meson or proton-nucleon interactions; the 
second type of special source is one for which the photo- 
meson (or proton-nucleon) optical depth is high. We be- 
gin by summarizing in Sec. O various speculations regard- 
ing the existence of such sources, for which we have no 
observational evidence to date. This section is particu- 
larly relevant for future neutrino observational programs. 
We then turn to a detailed discussion of the derviation 
of the bound and of its validity and robustness. 

In Sec. [II and Fig. nl we summarize the existing data 
on the cosmic energy spectrum, including the limited in- 
formation at very high energies that is derived from the 
Fly's Eye 111, AGASA |l|l, and Yakutsk H experi- 
ments. We show in this section that cosmic-ray observa- 
tions imply a a robust upper limit on neutrino fluxes in 
the energy range of ^ 10^^ eV to ^ 10^° eV. What is ex- 
pected at very high energies, > 10^" eV, and at relatively 
low energies, < 10^^ eV? We argue in Sec. IV that it is 



plausible to extrapolate the Waxman-Bahcall bound to 
neutrino energies that correspond to protons more ener- 
getic than the expected Greisen-Zatsepin-Kuzniin cutoff 
at 5 X 10^^ eV. In Sec. Q we summarize the observa- 
tional evidence for the conventional viewpoint that most 
cosmic rays below 10^* eV are heavier particles of Galac- 
tic origin. Thus, the proton flux that is to be used in 
deriving the upper limit to the neutrino flux is much less 
than the total observed cosmic ray flux at energies below 
10"'^^ eV. We discuss the important question of whether 
sources that are not included in the Waxman-Bahcall 
bound derivation may exist in this energy range. 

It has been common practice for some time to mo- 
tivate proposed observatories to search for high energy 
neutrinos by considering the predictions of AGN mod- 
els for neutrino production for photo-meson interactions 
that also produce the gamma-ray background. We show 
in Sec. VI that the large set of previously published AGN 
models for neutrino production exceed the Waxman- 
Bahcall bound by typically two orders of magnitude and 
therefore are not appropriate theoretical models upon 
which to base proposed observatories. We present es- 
timates for the upper limit neutrino event rates photo- 



meson and nucleon-meson interactions that may be ob- 
served in the AMANDA, ANTARES, Auger, ICECUBE, 
NESTOR, and OWL detectors. The previously pub- 
lished AGN models based upon photo-meson interactions 
are also inconsistent with gamma-ray observations that 
show rapid time variability. We also discuss in this sec- 
tion the question of whether some AGN's models may 
have a high photo-meson optical depth. 

Can one increase the rate at which high energy pro- 
tons are generated beyond the value contemplated in the 
Waxman-Bahcall bound? If so, one could obviously raise 
the experimentally important bound on the neutrino flux. 
The only way of exceeding the bound is by hypothesizing 
the existence of high energy proton accelerators which ex- 
ist in environments that do not allow protons to escape.[j 
Can one fine-tune the conditions in which luminous astro- 
nomical sources are imagined to exist so that neutrinos 
and photons escape, but protons do not leave the system? 
Strong magnetic fields obviously provide a potential way 
of confining protons, but not neutrons, and we discussed 
this possibility briefly in Ref. |1|]. In Sec. VII, we consider 
some specific scenarios for avoiding the Waxman-Bahcall 
limit by invoking strong mangetic fields and show by ex- 
ample that magnetic field confinement is not a plausible 
way of aoviding the limit. 



In Sec. VIII we show that it is a good approxima- 



tion for GRB's to assume that the high energy neutrino 
flux is dominated by photo-pion interactions that proceed 
through the A resonance. The contributions from non- 
resonant itneractions are shown to be small, comfirming 
our previous estimates P,H, according to which ^ 20% 
of the proton energy is converted in GRBs to pions. We 



summarize our results in Sec. IX 



II. EXCEED, NOT VIOLATE, THE WB BOUND 

Does the Waxman-Bahcall bound imply that neutrino 
fluxes from all conceivable astronomical sources must lie 
below this limit? Are all possible neutrino sources that 
will be searched for with very large area neutrino tele- 
scopes subject to this bound? 

As we discussed in the Introduction, there are two 
(speculative) ways to exceed the limit without violat- 
ing the constraint imposed by the observed cosmic ray 
background. The WB bound refers, see Sec. [II, only to 



sources that are optically thin to proton photo-meson and 



*We point out in Ref. |l| that one can imagine "neutrino 
only" sources that are optically thick to proton photo-meson 
interactions and from which protons cannot escape. Sources 
of this kind could in principle contribute a flux in excess of 
the Waxman-Bahcall bound, but there is, by construction, no 
observational evidence (from baryons or high energy photons) 
for their existence. 



proton-nucleon interactions and from which protons can 
escape. Thus, the first 'way out' which was also discussed 
in our original paper O], is to assume the existence of 
sources which are optically thick to proton photo-meson 
or proton-nucleon interactions. We referred in |l| to such 
sources as "neutrino only factories" or "hidden core mod- 
els" since they are not motivated by measurements of the 
cosmic-ray flux or by any electromagnetic observations. 

Examples of optically thick scenarios include neutrinos 
produced in the cores of AGNs (rather than in the jets) 
by photo-meson interactions |14[ , or via p — p collisions in 
a collapsing galactic nucleus |15| or in a cacooned black 
hole | 16|]. The optimistic predictions of the AGN core 
model]M| have already been ruled out by the AMANDA 
experiment [17| (see Fig. ^. 

The second type of source that may exceed the WB 
bound is one in which neutrinos are produced by some 
mechanism other than photo-meson or proton-nucleon 
interactions. We list below some of the hypothesized 
sources of this type. 

• Decays of supermassive, dark matter parti- 
cles |l8,n9[ might produce high energy neutrinos 
without violating the cosmic ray bound. 

• Topological defects (see, e.g., Refs. [po|"P5[) might 
produce particles which decay, among other chan- 
nels, to neutrinos with only a small associated cos- 
mic ray proton flux. 

• Superheavy relic neutrinos |2qj27|| which interact 
with the cosmic neutrino background have been 
proposed as one way of producing ultra-high en- 
ergy cosmic rays. 

• Ultrahigh-energy photons at large redshifts p^ ] 
could in principle produce a flux of ultrahigh- 
energy neutrinos. 



III. COSMIC RAY OBSERVATIONS AND THE 
WB BOUND 

In this section, we summarize the observations and cal- 
culations that lead directly to the upper bound on high 
energy neutrino fluxes. 

Figure ^ shows the cosmic ray fluxes measured by the 
Fly's Eye ^, AGASA |l|, and Yakutsk [|| experi- 
ments. The smooth curve shown in Fig. M was used by 
us [l| in setting a conservative upper bound on the high 
energy neutrino fluxes. We now explain why the upper 
bound is robust and conservative. 

The smooth curve was computed assuming that in the 
nearby universe (z = 0) the energy production rate is 
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FIG. 1. The observed high energy cosmic ray flux. Mea- 
surements are shown from the Fly's Eye Q, AGASA Q, 
and Yakutsk |l3] detectors. The smooth curve, computed 
from Eq. (nl), was used by Waxman and Bahcall |1| to com- 
pute the upper bound on high energy astrophysical neutrino 
sources from p — 7 interactions. 

It is possible that the energy generation rate given lo- 
cally by Eq. (|l|) increases with redshift. In particular, 
we have explained in Ref. |l| how cosmic rays observed 
at earth to have energies in excess of 10^^ eV must have 
originated at small redshifts because of the large energy 
loss rate at these high energies. In order to establish a 
conservative upper limit, we assumed that the local rate 
given in Eq. (Q) evolves with redshift at the maximum 
rate observed for any astronomical population, i.e., the 
evolutionary rate exhibited by the quasars p9|-|3l|. We 
also included the adiabatic energy loss due to the expan- 
sion of the universe. 

Figure |l| shows that the smooth curve which we have 
used to estimate the cosmic ray flux above 10^^ eV is a 
conservative (i.e., high) estimate of the observed rate.| 

What is the neutrino bound that results from the ob- 
served cosmic ray flux? Figure g shows the numerical 
limit that is implied by the cosmic ray observations. The 
upper horizontal curve is computed by assuming that the 
cosmic ray sources evolve as rapidly as the most rapidly 
evolving known astronomical sources. This very conser- 
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^ We note that Fig. hi shows that the highest energy point 
measured by the AGASA experiment could be interpreted to 
suggest (with ~ Icr significance) that the cosmic ray gener- 
ation rate aX E > 10^" eV is twice the rate obtained from 
our smooth curve generated by Eq. (|l|), implying that the 
upper bound might be underestimated by a factor of two at 
~ 10^^ eV. However, the higher rate of generation is not ob- 
served by the Fly's Eye and Yakutsk experiments, and even 
if correct would imply only a small correction to the upper 
bound at this energy. 



vative limit is what we shall mean when in the follow- 
ing we refer to the 'Waxman-Bahcall bound.' The lower 
curve is computed assuming that the number density of 
cosmic ray sources at large distances is the same as in 
the local universe. We will discuss the implications of the 



bound for previously published AGN models in Sec. Vl 
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FIG. 2. The Waxman-Bahcall (WB) upper bound on muon 
neutrino intensities (j^^ -|- iv^). The numerical value of the 
bound assumes that 100% of the energy of protons is lost 
to tt"*' and tt" and that the tt^ all decay to muons that also 
produce neutrinos. The WB upper bound exceeds the most 
likely neutrino flux by a factor of 5/r for small optical depths 
r. The upper solid line gives the upper bound corrected 
for neutrino energy loss due to redshift and for the maxi- 
mum known redshift evolution (QSO evolution, see text). In 
what follows, we will refer to this conservative upper curve 
as the "Waxman-Bahcall bound." The lower solid line is ob- 
tained assuming no evolution. The dotted curve is the maxi- 
mum contribution due to possible extra-galactic component of 
lower-energy, < 10^^ eV, protons as first discussed in ll| (see 
Sec. Mfor details). The dash-dot curve shows the experimen- 
tal upper bound on diffuse neutrino flux recently established 
by the AMANDA experiment [[L7|. The dashed curves show 
the predictions of the GRB fireball model [@,[y,K2|. 



IV. BEYOND 10^" EV 

The theoretical curve in Fig. |l| shows the predicted 
decrease above 5 x 10^^ eV in the observed cosmic ray 
flux due to interactions with the cosmic microwave back- 
ground [Q. The existing obervations do not clearly 
demonstrate this decrease and therefore many authors 
have speculated that there may be a new source of ultra- 
high energy neutrinos beyond 10^° eV. 

In our original discussion (see Fig. 1 of Ref. ||l| or Fig. g 
of the present paper), we did not extend the upper bound 
to energies beyond 10^° eV because we judged that too 
little is presently known about the cosmic ray flux in that 
region. However, as a useful guide to designing experi- 



ments, it is a reasonable expectation that the production 
rate will not exceed by a large factor the flux predicted by 
an extension of the E~'^ production rate used in setting 
the limit at observed proton energies of 10^^ eV. 

In our view, it is reasonable to extend the Waxman- 
Bahcall limit beyond 10^° eV by simply extrapolating the 
horizontal line in Fig. 0to higher energies. The validity of 
this extrapolation will be tested by future measurements 
of the spectrum of ultra-high energy cosmic rays, mea- 
surements that will determine observationally whether 
the spectrum exhibits the predicted Greisen-Zatsepin- 
Kuzmin (GZK) cutoff at high energies. 

We have suggested elsewhere Q that the clustering 
of sources of ultra-high energy (> 10^°) cosmic rays may 
give rise to more events beyond the GZK cutoff than 
expected on the basis of a homogeneous distribution of 
sources. Other authors (see, for example, the references 
cited in Sec. O) have suggested exotic sources of ultra- 
high energy cosmic rays that are not described by the 
E^'^ proton generation rate. If any of these scenarios is 
correct, then the upper bound shown in Fig. g will have 
to be raised above the simple extrapolation we currently 
advocate as a useful guide. However, the existing cosmic 
ray data suggest that the factor by which the WB upper 
bound is conservative, 5/t, will more than compensate 
for any increase in cosmic ray flux beyond the GZK limit. 



V. LOWER ENERGY: E < 10^** EV 

Figure ^ shows that our smooth curve for the extra- 
galactic proton cosmic ray flux falls below the total ob- 
served flux for energies less than 10^^ eV. We have as- 
sumed (see Eq. |l|) that the generation rate of extragalac- 
tic protons is proportional to E~'^, where E is the proton 
energy. This E^'^ dependence is produced generically by 
the Fermi mechanism for accelerating high energy cosmic 
rays in shocks pq| . If the lower energy cosmic rays were 
protons from extragalactic sources, then (as we noted ear- 
lier M) one could raise the upper bound for astrophysical 
neutrinos at these lower energies. 

Could the lower energy cosmic rays be primarily ex- 
tragalactic protons? Unfortunately, the answer is nq^. 
The available observational evidence suggests that only 
a small fraction of the cosmic-ray flux in the energy range 
of 10^^ eV to 10^^ eV is composed of protons. Direct (bal- 
loon) composition measurements at 10^"* eV ||3q| show 



"'■In the original version of their paper, astro-ph/9812398 
vl, Mannheim, Protheroe, and Rachen assumed that 100% 
of the lower energy cosmic rays could be extragalactic pro- 
tons. In the third version of their paper, v3, they acce pted 
the argumen t presented in the first version of our paper |iep-| 



ph/9902383 vl, and repeated here, that only a small fraction 
of these lower energy cosmic rays are extragalactic protons. 



that the fraction of cosmic-ray flux composed of pro- 
tons at this energy is ~ 20%. Air-shower and cosmic-ray 
tracking detectors measurements indicate gj that the 
proton fraction decreases in the energy range of 10^^ eV 
to 10^6 eV. In fact, the Fly's Eye and AGASA experi- 
ments support |38|] a composition strongly dominated by 
(consistent with 100%) heavy nuclei at 10"'^'^ eV. More- 
over, both the Fly's Eye @ and the AGASA |§ collab- 
orations have reported, for energies less than 3 x 10^^ eV, 
a small but statistically significant enhancement of the 
cosmic ray flux near the Galactic plane. This enhance- 
ment is expected for a Galactic, but not an extragalactic, 
origin for the cosmic rays in this energy domain. 

In summary, the observational evidence is that most 
of the observed cosmic rays in the energy range 10^"' eV 
to 10^^ eV are not protons and therefore the total flux of 
cosmic rays cannot be used to raise the upper bound for 
neutrinos in the energy range 10^*^ eV to 10^^ eV. The 
NASA satclhte ACCESS is expected to provide accurate 
measurements of the composition of the cosmic rays at 
energies up to ~ 10^^ eV sometime in the next decade. 

Assuming, conservatively, that ^ 10% of the cosmic 
rays in this energy region are protons, the neutrino bound 
may be raised at energies E^, < 10^^ eV (Since the cosmo- 
logical model we have used, see Eq. |l| and Fig. m, accounts 
for ^ 10 % of the total cosmic ray flux at 10^ eV, and a 
progressively larger fraction at higher energies, one can- 
not increase the upper bound on neutrino fluxes at ener- 
gies > 10"'^^ eV). Since the energy density in cosmic rays is 
approximately proportional to E~^ for energies less than 
10^* eV, the actual extragalactic neutrino flux could ex- 
ceed the most stringent Waxman-Bahcall limit by a fac- 
tor of (10^^ e\[ /By) for neutrino energies Ey < 10"^^ eV. 

Figure g illustrates with a curved dotted line the 
maximum contribution from unrecognized extragalactic 
sources of protons, i.e., sources of extragalactic protons 
that do not contribute significantly to the observed cos- 
mic ray flux at 10^^ eV. The 'unknown source' contribu- 
tion is obtained by multipling the proton generation rate 
that is proportional to E~'^ by an analytic approxima- 
tion to the maximum additional contribution allowed at 
lower energies by the existing experiments. Thus 

max. proton rate oc E^'^ [l + Q.l{{lO'^'^eV)/E)] , (2) 

where E is the proton energy. The recent experimental 
bound established by the AMANDA experiment already 
rules out a contribution from unrecognized extragalac- 
tic sources of protons at the maximum level allowed by 
present cosmic-ray observations. 

We note that including a contribution to the neu- 
trino flux due to neutrino production by photo-meson 
interactions on heavy nuclei will not affect the neutrino 
bound. For heavy elements, the cross section for photo- 
dissociation is higher than, and its energy threshold is 
lower than, for photo-meson production. Thus, heavy el- 
ements will dissociate before losing a significant fraction 
of their energy to neutrino production. 



VI. AGN MODELS 

What predictions have been made for the neutrino 
fluxes from AGN jets? Are the fluxes large enough to 
be measurable in a practical detector? 

There are two classes of models that have been consid- 
ered as explanations for the emission of radiation from 
AGN jets, ones in which radiation is due to electromag- 
netic processes [Q and others in which emission of radi- 
ation involves the acceleration of very high energy pro- 
tons [[l2|-p4[- For the first class of models, no significant 
flux of high energy neutrinos is expected. The second 
class of models has been normalized by postulating that 
the jets produce the observed gamma-ray background via 
photo-meson processes on high energy protons, which 
produce tt^'s that decay into gamma rays. The photo- 
meson process, if it produces the gamma-ray background, 
would also produce a large flux of high energy neutrinos 
via charged meson decay. Therefore, the second class of 
models, involving photo-meson interactions on high en- 
ergy protons, has received a lot of attention, especially 
from particle experimentalists. 

The previously published AGN models that are dis- 
cussed in the following subsection have been used to 
estimate the expected event rates in the AMANDA, 
ANTARES, Auger, ICECUBE, NESTOR, and OWL de- 
tectors p-10|. However, all of these models exceed by 
more than an order of magnitude the maximum neutrino 
flux permitted by the WB bound (see Fig. S). 



A. Previously published models 

Figure compares the Waxman-Bahcall upper bound 
to predictions by various proton acceleration models for 
AGN jets. The previously published models illustrated 
in Fig. ^ are inconsistent with the cosmic ray limit; they 
typically predict fluxes one to two orders of magnitude 
above the Waxman-Bahcall upper bound. 

The expected detection rate in a km^ detector (e.g., 
ANTARES, ICECUBE, or NESTOR) is less than one 
event per year from AGNs with spectral neutrino energy 
shapes like the models P97, IIZ97, and M95B shown in 
Fig. y if a neutrino flux consistent with the Waxman- 
Bahcall bound is assumed. This rather pessimistic result 
can be obtained from Tables III- VI of Ref. [^ by di- 
viding by 30 the detection rate given in those tables for 
the model denoted by AGN-M95. The same model is la- 
beled 'M95B' in Fig. H, where one can see that the flux 
predicted from M95B is about a factor of 30 above the 
cosmic ray upper bound. 

The neutrino event rate in the Auger detector' pjlCJ] 
is less than or of order 1 event per year above a neutrino 
energy of 10^^ eV for optically thin sources, i.e., sources 
that satisfy the WB bound. The rate may be two orders 
of magnitude larger if the OWL/AIRWATCH detector is 
built according to the preliminary specifications M . 
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FIG. 3. The Waxman-Bahcall (WB) upper bound on muon 
neutrino intensities (i^^j + v^) (solid line) compared to pre- 
dictions of representative AGN jet models, taken from the 
earlier papers of Mannheim |44| (marked M95B in the fig- 
ure), Protheroe Q (P97), and Halzen and Zas jii) (HZ97). 
The AGN models were normalized so that the calculated 
gamma-ray flux from tt" decay fits the observed gamma-ray 
background. The AGN hidden-core conjecture (S92), to 
which the WB upper bound does not apply due to high 
photo-meson optical depth of the source, is taken from 
[ [l4[ . Note, that this conjectrure is already ruled out by the 
AMANDA upper bound Q. 

Since p — 7 reactions give rise to high energy gamma 
rays (via tt'^ decay) as well as to high energy neutrinos 
(via charged pion decay), the flux of neutrinos is pro- 
portional to the flux of gamma rays. Since the models 
with high neutrino fluxes that are shown in Fig. g were 
all chosen by their authors so that they could account 
for the observed gamma ray flux, it is clear that the still 
allowed AGN models (with the 30 times lower neutrino 
flux) can only account for a few percent of the gamma 
ray background. 



B. Are AGNs optically thick to photo-meson 
interactions? 

The WB bound applies only to sources which are op- 
tically thin to photo-meson interactions^. As mentioned 
above, previously published AGNjet models satisfy this 
requirement. One may argue Bq], however, that new 



^The photo-meson optical depth of high-energy protons in 
GRBs is discussed in |llp|. It is shown there that the obser- 
vation of high energy, ~ 100 MeV, photons in GRBs implies 
that the sources are not optically thick to photo-meson in- 
teractions. In the context of Fireball models for GRBs, the 
optical depth to nucleon-meson interactions is also small 



models may be constructed to have large optical depth 
and nevertheless fit the available data. 

We discuss this possibility below. We first explain why 
published AGN models predict small optical depth, and 
then discuss whether new, high optical depth models can 
be constractcd. 

The optical depth to photo-messon interactions for 
protons of energy Cp can be related to the optical depth 
for pair production of photons of much lower energy, e^. 
The threshold relation for photon-meson interactions is 
tp e^ « TTT.^ nip and the threshold relation for pair pro- 
duction is e' e-y — 2TOg. Thus, as first pointed out in [nl, 
a photon of energy e.y that causes a photon-meson inter- 
action with a proton of energy tp could also pair-produce 
with a photon of energy [2mg/(m7rm.p)]ep = 4 x 10~^ep. 
Taking account of the ratio between photo-meson and 
pair production cross sections, one finds that 



'tpV^P 



i^p) 



5 X 10" 



T^^{ej ^ 4 X 10 



ep)- 



(3) 



Observed AGN photon energy distributions typically fol- 
low a power-law, dn^/de^ oc e~^. For such photon spec- 
tra, one can easily show |l|] that t-^^ cx e^, and hence 
that 



,p(6p) «2r,,(6, = 10GeV)(^^) 



(4) 



Emission of ~ 1 TeV photons from "blazars," AGN 
jets nearly aligned with our line of sight, is now well 
established p3, and there is evidence that the high- 
energy photon spectrum extends as a power-law at least 
to ~ 10 TeV p^. The high energy photon emission 
is the main argument used pfl] in support of the hy- 
pothesis that high-energy emission from blazars is due 
to pion decay rather than inverse Compton scattering. 
However, the observation of > 1 TeV photons implies, as 
shown by Eq. (|3|), that the jet optical depth to photo- 
meson interaction is very small. Indeed, one can read- 
ily see from Fig. 4 of Ref. pM that high-energy photon 
data of northern hemisphere blazars requires small pair- 
production optical depth at photon energy < 0.1 TeV. All 
protons blazar models shown in the figure have r^^ < 1 
at this energy, i.e. r^p < 1 for Sp < 10^^ eV. 

TeV emission is observed from the nearest blazars, 
which are hence relatively low-luminosity blazars. It 
had recently been argued by Mannheim, Protheore and 
Rachen (MPR, |4^), that the emission of high energy, 
~ 1 TeV, photons may be suppressed in high luminosity 
blazars, for which models with high optical depth may 
therefore be constructed. The most intense high-energy 
gamma-ray source is 3C 279. Figure 6 of Ref. |4^], which 
is reproduced here as Fig. 0, shows that the measured 
optical depth for pair production is at most r^..^ ~ 1 at 
e^ = 10 GeV. Combining this observation with Eq. (0) 
we see that r-yp(ep) < 1 for e^ < 10^® eV. Making the ad 
hoc assumption, that the pair-production optical depth 
for high luminosity blazars like 3C 279 exceeds unity for 
photons above 10 GeV, i.e. just above the highest energy 



for which data are available (see Fig. ^) , MPR concluded 
that the AGN neutrino flux may slightly exceed the WB 
bound at ~ 10^^ eV (see dotted curve in their Fig. 5a). 
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FIG. 4. The emission spectrum of 3C 279 over 15 decades 
in frequency, adatped from Ref . |47| . The measured fluxes are 
plotted as small squares in the figure. The flare in June 1991 
was studied by many difi^erent groups observing in a large 
number of wavelength domains. Out to the last measured 
point at about 10 GeV, the emission spectrum shows no evi- 
dence for turning over due to a large optical depth. The solid 
curves are the predictions for emission from a uniform rela- 
tivistic moving sphere; more complex models are discussed in 

In order to avoid the constraint given in Eq. (0), and 
construct AGN models with significant optical depth to 
photo-meson interaction, one must argue that the spec- 
tral distribution of the background photons (for photo- 
meson and pair-production intercations) deviates from 
a dn-f/dcj ex e~^ power-law. Since the observed spec- 
tral distribution does follow an e;^^ power-law, one must 
argue that the background radiation with which high en- 
ergy photons and protons may interact is not the ra- 
diation associated with the observed photon flux. This 
may indeed be the case for jets which expand with high 
Lorenz factor F. In this case, an ambinet isotropic ra- 
diation field, Uamb.i through which the jet expands, will 
dominate pair-production and photo-meson interactions 
as long as its energy density exceedss a fraction F^** of the 
energy density associated with the observed radiation, 
Uobs., which is produced within the jet. The presence 
of such radiation field can not be ruled out. However, 
the energy density of such ambinet radiation is limited 



to values not much exceeding Uamb./uo. 



bs. 



T- 



for Uamb./uobs. > T"^ invcrsc Compton scattering emis- 
sion of jet electrons would dominate over the synchrotron 
emission which produces the observed radiation in the ra- 
dio to UV bands. 

In order to obtain the maximum possible optical depth, 
MPR have therefore adopted Q] the following ad hoc 
assumptions: (1) The pair production optical depth in 
luminous blazars, e.g. 3C 279, exceeds unity for photons 
above 10 GeV; (2) The optical depth is due to ambient 
radiation (which is not detected); (3) The spectrum of 
the ambinet radiation deviates from a e~^ power law in 
such a way that its energy density increases by an order 
of magnitude at photon energies just below those cor- 
responding to the pair-production threshold of 10 GeV 
photons. Adopting these implausible assumptions, MPR 
concluded that the AGN neutrino flux may exceed the 
WB bound by a factor ^ 2 at ^ lO^'* eV (see thick dahsed 
curve in their Fig. 5a). 



VII. CAN ONE EVADE THE UPPER BOUND 
WITH THE AID OF MAGNETIC FIELDS? 

In Sec. Ill of Ref. Q we considered various ways that 
one might try to avoid the upper limit on the flux of high 
energy neutrinos by invoking magnetic flelds. In particu- 
lar, we discussed the possibility that the cosmic-ray den- 
sity observed at Earth is lower than the average universe 
cosmic-ray density, due to confinement of protons by ma- 
gentic fields surrounding the cosmic-ray sources, or due 
to large-scale structure fields. We have shown that, given 
observational constraints on such magentic fields, mag- 
netic proton confinement can not affect the high-energy 
cosmic-ray fiux, and therefore can not affect the neutrino 
bound. 

In Sec. IIIc of Ref. [Q we have discussed in de- 
tail the possible effects of large-scale structure magnetic 
fields, and demonstrated that even if such fields are 
built to equipartition levels, protons of energy exceed- 
ing ~ 10^^'^ eV can not be confined to large-scale struc- 
tures over a Hubble time. Large-scale structure magnetic 
fields can not therefore affect the WB bound at energy 
> lO^'^ cV. It should be pointed out here, that had pro- 
ton confinement by large scale structure magnetic fields 
been possible, which might be the case for low energy 
protons, it would have most likely implied that the neu- 
trino upper bound should be lower, rather than higher at 
Ejy ^ 10^^ eV. This is due to the fact that confining high- 
energy cosmic-rays to high density regions would imply 
that the local cosmic-ray density is higher than the uni- 
verse average. Since, on the other hand, neutrinos can 
not be confined and their density is uniform, their flux 
would be bound to a lower value than implied by our 
assumption of a uniform cosmic-ray density. 

Recently, it has been argued by Mannheim, Protheore 
and Rachen (MPR, El]), that while high energy pro- 
tons can not be confined to the vicinity of their sources 



by magnetic fields, adiabatic losses of protons escaping 
magnetic halos surrounding AGN jets may lead to signif- 
icant decrease in proton energy for ep < 10^^ eV, which 
may lead to modification of the bound for €^ < 10^^ eV. 
Clearly, even if this claim is correct, proton AGN mod- 
els are still in contradiction with the Waxman-Bahcall 
bound, since the bound is not affected at ei, > 10^^ eV 
(see Fig. H). Moreover, the key arguments made by 
MPR in deriving the claimed adiabatic energy loss are 
flawed. First, the magnetic field model adpoted by MPR 
is in conflict with observations. While the jet mag- 
netic fleld inferred from observtaions is high, tens of 
/iG (e.g. |52]), depolarization measurements |53] imply 
small halo magnetic field, 0.1-1 //G at the central halo 
region, r ~ 10 kpc, and much smaller at large distances, 
r r^ 100 kpc. MPR assume a spherically symmetric field 
structure with amplitude decreasing as r~^ from 30 /uG 
r = 10 kpc to 0.3 fj,G at r = 100 kpc. Second, MPR as- 
sume that high energy protons are produced in the cen- 
tral region of the jet, where the magnetic field is high, and 
then are confined to jet plasma which expands adiabati- 
cally to the halo plasma conditions of low magnetic field. 
This is an implausible description of the plasma evolu- 
tion and of proton confinement. Acceleration of protons 
to high energy is typically expected to occur at the outer 
edge of the jet at the strong shock. Moreover, since the 
jet is narrow, typical opening angles are 10 degrees, the 
typical jet structures are ~ 1 kpc in scale, and neutrons 
produced by photo-meson interactions leading to i^ pro- 
ductiuon escape the jet before decaying. 



VIII. CONTRIBUTIONS FROM 
NON-RESONANT INTERACTIONS 

The high energy-neutrino flux from GRBs was cal- 
culated by Waxman & Bahcall [g|j^] using the "A- 
approximation," i.e. assuming that photo-meson produc- 
tion is dominated by interaction of protons with photons 
of energy close to that corresponding to the A-resonance. 
It has recently been argued by Miicke et al. [Q that for 
characteristic GRB photon spectrum, there is an addi- 
tional contribution to photo-meson production from in- 
teraction of protons with photons of energy much higher 
than that corresponding to the resonance, leading to sig- 
nificant deviation from results derived based on the A- 
approximation. We point out here, that the contribu- 
tion to neutrino production of non-resonant interactions 
is small, leading to only negligible modification of results 
obtained using the A-approximation. 

The GRB photon spectrum is weU fitted in the BATSE 
detector range, 30 keV-3 MeV, by a combination of two 
power-laws, dn^/de^ oc e~^ with different values of /3 at 
low and high energy [p5| . The break energy (where /3 
changes) in the observer frame is typically e°^' ~ IMeV, 
with /? ~ 1 at energies below the break and /3 ~ 2 above 
the break (the plasma emitting the GRB radiation ex- 
pands with Lorentz factor F ~ 300, and particle ener- 



gies in the plasma frame are smaller by a factor F~^ 
compared to observed particle energies). For protons of 
observed energy ep < 10^^ eV, the (observed) photon 
threshold energy for photo-meson interaction is higher 
than the break energy e°^' ~ IMeV B. Such protons 
therefore interact only with the steep part of the pho- 
ton spectrum, dn^/de^ ex e~^. Thus, for ep < 10"'^^ eV 
the contribution from interaction with high energy pho- 
tons, i.e. photons of energy well above the threshold and 
therefore well above the energy corresponding to the A- 
resonance, is negligible simply because there are very few 
photons of such high energy. 

Protons of energy ep ^ 10^^ eV may interact with 
photons well below the break, where the photon spec- 
trum is flatter, dn^/de^ oc e"^. In this case there may be 
significant contribution to photo-meson production from 
interaction of protons with photons of energy well above 
the threshold, and therefore well above the energy corre- 
sponding to the A-resonance, since the photon spectrum 
is flat over significant energy range, from the threshold 
energy to e°^' ~ IMeV. However, contribution from 
non-resonant interaction becomes significant only at very 
high energy: for protons of energy ep > 10^® eV photo- 
meson production of pions is increased by only a factor 
of '^ 2 compared to the production rate based on the A- 
approximation. Moreover, this increase has no effect on 
the neutrino production discussed by Waxman & Bahcall 
iy,y , since the high energy pions produced in interaction 
of protons with ep > 10"'^* eV lose their energy through 
synchrotron emission before decaying, and therefore do 
not contribute to the neutrino flux, which is strongly sup- 
pressed at observed neutrino energy > 10^^ eV P,|56|,p| . 



IX. DISCUSSION 

In this paper, we have shown that the Waxman-Bahcall 
upper limit on neutrino fluxes is robust for neutrino ener- 
gies above 10^^ eV. Cosmic ray experiments set a firm up- 
per limit on the flux of extragalactic protons, and hence 
on the flux of extragalactic neutrinos, that are produced 
by P + 7 OT p — p interactions in astronomical sources that 
are optically thin to photo-meson and nucleon-meson in- 
teractions on high energy protons. The upper solid curve 
in Fig. H, which is labeled WB bound, is a conservative 
representation of the limit imposed by the cosmic ray ob- 
servations and includes the maximum plausible redshift 
evolution of the sources of high-energy cosmic rays. The 
WB bound exceeds the most likely neutrino flux from 
p + 7 OT p — p interactions by a factor of 5/r, for small 
optical depths, r, to photo-meson and nuclcon-mcson in- 
teractions. 

The Waxman-Bahcall upper bound is about two orders 
of magnitude less than published predictions |4^-Q of 
neutrino fluxes expected from models for AGN jets that 
explain the gamma-ray background by photo-meson in- 
teractions on high energy protons (see Fig. 0) . Given the 



neutrino upper limit, photo-mesonic interactions in opti- 
cally thin sources like AGN jets can at most account for a 
few percent of the flux of the gamma-ray background [Q . 

The Waxman-Bahcall bound implies an upper limit 
to the expected event rate from optically thin sources 
that may be observed in the AMANDA, ANTARES, 
Auger, ICECUBE, NESTOR, and OWL/AIRWATCH 
detectors |-0. In km^ detectors like ANTARES, ICE- 
CUBE, or NESTOR, one expects to detect less than or of 
order 1 neutrino event per year from optically thin AGNs 
with neutrino spectral energy shapes like the AGN mod- 
els P97, HZ97, and M95B shown in Fig. |. For the Auger 
detector, one expects less than or of order 1 neutrino 
event per year above 10^^ eV from all sources satisfying 
the WB bound. The event rate could be ^ two orders 
of magnitude larger for the OWL detector if it performs 
according to the preliminary specifications. These esti- 
mates are upper limits and do not include the factor of 
5/t by which, for small optical depths r, the WB bound 
exceeds the plausible expected neutrino flux. The rate 
estimates given in this paragraph are much less than the 
estimates given in the published papers describing the ca- 
pabilities of the neutrinos detectors, since in those papers 
AGN models like those shown in Fig. pi were assumed in 
the calculations. 

At neutrino energies below 10^^ eV, the flux of extra- 
galactic neutrinos from optically thin sources may reach 
as high as the dotted curve in Fig. g. In order for this 
curve, labeled 'max. extra-galactic p' in Fig. 0. to ap- 
ply, a large fraction of the cosmic rays with energies be- 
low 10^^ eV must be protons of extra-galactic origin (cf. 
Eq. H). For cosmic rays more energetic than 10^° eV, 
there is not yet an accurate determination of the cosmic 
ray flux. Hence there is no way of rigorously extending 
the WB bound for neutrino energies > 10^^ eV. We pro- 
pose that horizontal line in Fig. simply be extrapolated 
to higher energies as a plausible bound if the source spec- 
trum continues to fall off at ultrahigh-energies as E~^. 

The large area neutrino experiments y-|7| that are cur- 
renty being constructed are designed to detect neutrinos 
with energies < 10^^ eV, corresponding to protons with 
energies less than or of the order of 10^^ eV. The WB 
limit is constructed by normalizing to the observed flux 
of cosmic ray protons at 10^^ eV and extrapolating to 
lower energies using the E~^ source spectrum. The ex- 
trapolation is plausible but not rigorous for the neutrino 
energy domain of the ANTARES, BAIKAL, ICECUBE, 
and NESTOR experiments. 

There are two hypothetical ways of exceeding the 
Waxman-Bahcall limit without violating the constraint 
imposed by cosmic ray observations. First, the neutrinos 
could be produced in sources that are optically thick to 
photo-nucleon or nucleon-nucleon interactions. Second, 
the neutrinos could be produced by processes that do not 
give rise to high energy cosmic rays, such as the decay 
of dark matter particles, topological defects, superheavy 
relic neutrinos, or ultrahigh-energy photons. So far, there 
is no obervational evidence supporting any of these pos- 



siblities. 

We have investigated a number of ad hoc scenarios 
invented to try to find ways of violating the Waxman- 
Bahcall bound. None of the suggested scenarios, in- 
cluding various ideas involving magnetic fields, provide a 
physically self-consistent mechamism for raising the up- 
per limit to the neutrino flux implied by a straightforward 
interpretation of the cosmic ray observations. 
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